How accurate is Poisson-Boltzmann theory for monovalent ions near highly charged 

interfaces? 
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Surface sensitive synchrotron x-ray scattering studies were performed to obtain the distribution 
of monovalent ions next to a highly charged interface. A lipid phosphate (dihexadecyl hydrogen- 
phosphate) was spread as a monolayer at the air-water interface to control surface charge density. 
Using anomalous reflectivity off and at the I/3 Cs"*" resonance, we provide spatial counterion (Cs+) 
distributions next to the negatively charged interfaces. Five decades in bulk concentrations are inves- 
tigated, demonstrating that the interfacial distribution is strongly dependent on bulk concentration. 
We show that this is due to the strong binding constant of hydronium H3O''' to the phosphate 
group, leading to proton transfer back to the phosphate group and to a reduced surface charge. The 
increase of Cs"*" concentration modifies the contact value potential, thereby causing proton release. 
This process effectively modifies surface charge density and enables exploration of ion distributions 
as a function of effective surface charge-density. The experimentally obtained ion distributions are 
compared to distributions calculated by Poisson-Boltzmann theory accounting for the variation of 
surface charge density due to proton release and binding. We also discuss the accuracy of our 
experimental results in discriminating among possible deviations from Poisson-Boltzmann theory. 
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I. INTRODUCTION 



The theoretical determination of ion distributions in 
aqueous solutions was initiated almost a century ago by 
Guoy 1] and Chapman who applied the Poisson- 
Boltzmann (PB) theory to calculate the spatial distri- 
bution of monovalent ions near a uniformly charged in- 
terface. Ever since their seminal work, the topic remains 
central in statistical mechanics, physical chemistry, and 
biophysics 3] . The original PB theory is a mean field the- 
ory with some simplified assumptions such as, point like 
particles and uniform surface charge density. To account 
for the finite ionic radius. Stern introduced phenomeno- 
logically a layer of ions at the charged interface with a dif- 
ferent dielectric constant, the Stern layer 1^. The effect 
of excess salt concentration and the resulting screening 
was extended by Debye-Hiickel . Grahame generalized 
the Gouy-Chapman theory to multivalent ions [g ■ Subse- 
quently, more refined theories and numerical simulations 
were developed to incorporate short-range interactions, 
image charges, finite size ionic radius, and ion-ion cor- 
relations 0) 0, Q • More recently, modifications of 
PB theory have been developed to incorporate hydration 
forces 0, . Some first-principles calculations of 
surface-tension for amphiphilic monolayers assume a PB 
theory with one or more layers of varying dielectric con- 
stant 14]. 

Experimental support for the validity of PB theory 
was provided by electrokinetic, visco-electric effects, and 
other techniques (see Ref. flSj for a review). McLaughlin 
and collaborators have shown good agreement be- 
tween C-potentials computed from PB theories and elec- 
trophoretic measurements in lipid vesicles. Other tech- 
niques, such as radiotracer experiments [tRIT^Ii^ . x-ray 
reflectivity polf2lll2^ . or infrared spectroscopy j23l| allow 



the determination of the total amount of ions in the im- 
mediate vicinity of a charged interface. It is noteworthy 
that all the experimental data for monovalent ions (at 
moderate salt concentrations < 0.1 M) outlined above 
are adequately described by the Guoy-Ghapman theory 
(with the generalization of excess salt) with no need for 
further corrections [2lll2^ . A close inspection, however, 
shows the agreement between theory and experiment is 
either based on fitting variables such as surface charge 
or interfacial dielectric values, that are not known in 
advance or/and based on integrated quantities. As an 
example, it has been shown recently that the degree of 
proton dissociation of arachidic acid spread on a sodium 
salt solution is adequately described by PB theory, but 
this agreement only involves the integral (over the en- 
tire space) of the sodium distribution Thus, 
local deviations that preserve the integral of the distri- 
bution (i.e., total number of ions) are not discriminated 
by these experiments. On the other hand, force measure- 
ments between two charged membranes separated by salt 
solutions, although well described by PB theory at large 
distances, show strong deviations at short distances (1-2 
nm) 0, [iJl . The origin of these hydration forces is still 
controversial. In some cases, it has been suggested to ex- 
tend PB theory to incorporate the restructuring of wa- 
ter, resulting in ion distributions that deviate from PB at 
short distances from the interface [T^l2^l26| . It is there- 
fore imperative to determine the ion distribution itself to 
establish the degree of accuracy of PB theory. 

In this manuscript, we experimentally determine the 
distribution of Cs+ ions next to charged interfaces by 
anomalous x-ray reflectivity |23, ll^ for several decades 
in bulk cesium concentrations. We compare our results 
with the predictions of PB theory and present a dis- 
cussion on the sensitivity of the experimental results to 
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quantify the magnitude of possible deviations from the 
PB distribution. This kind of investigation has become 
feasible only with the advent of the second generation x- 
ray synchrotron sources with novel insertion devices (i.e. 
undulator) and improved optics, which readily produce 
variable-energy x-ray beams with brilliances capable of 
detecting a single atomic-layer even if not closely packed. 

The manuscript is organized as follows: In Section ITTI 
we review several theoretical results relevant to this 
study. In Section IIIII the experimental details are pro- 
vided. In Section IIVI the x-ray reflectivity (XR) and 
grazing incidence x-ray diffraction (GIXD) measurements 
are presented with detailed structural analysis. The ex- 
perimental results are compared with the theoretical PB 
theory in Section The implications of the present 
study are discussed in Section Ivll 



II. PB THEORY INCLUDING PROTON 
TRANSFER AND RELEASE 

We briefly review the PB theory of an aqueous solu- 
tion containing monovalent ions (1:1 positive and nega- 
tive electrolytes) of average bulk concentration rif, in the 
presence of an ideally flat charged interface of known sur- 
face density (Tq = —s-uq, a.t z = Q. As the surface charge 
is uniformly distributed, the electric potential ip depends 
only on the distance from the interface z, and the Boltz- 
mann distribution is 

p{z) = enb[e^'='^(^)/'=«'^ - e'='^(^)/'=«^] = -2en6 sinh (/.(z), 

(1) 

where fc^ is the Boltzmann constant and 4>{z) = 
eip{z)/kBT . The ion distribution is calculated self- 
consistently from the Poisson equation 
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where \d = (eA:BT/87re^n;,)^/^ is the Debye screening 
length. Equation |(2J) can be solved analytically, yielding 
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FIG. 1: (A) Calculated monovalent ion distributions n^{z) 
near a negatively charged surface as obtained from PB the- 
ory, Eq. ||1J for different bulk concentration. Surface charge 
density ctq is one electron charge per 40 A"^, e — 80, and 
T — 292 K. Note that for this surface charge density the 
value of the distribution at the 2 = is practically a con- 

stant, n^(0) ~ 1^1% i^) The convolutions of the distribu- 
tions n^{z) using Eq. Q assuming two F values as indicated. 
For this strongly charged surface the convoluted distributions 
for different bulk ionic concentrations are practically indistin- 
guishable. 



where 7 = tanh[e-0(O)/4fcBr]=-AGc/An + 
{{Xgc/Xd? + 1)'/', and Agc = kBTe/2naoe is 
the Gouy-Chapman length. Figure ^A) shows cal- 
culated counter-ion distributions using Eq. Q for a 
fixed surface charge density, typical of fully deproto- 
nated closed packed phospholipids cto = — e/40 A^ for 
several salt concentrations. It is worth noting that 
for this surface density (Agc ~ 0.9 A) and bulk salt 
concentrations (Ad = 963.5 A and 9.6 A at 10~^ and 
0.1 M, respectively) << 1, and 7 « 1 - the 



concentration of ions next to the interface is independent 
of bulk concentration 

eksT z 
27re^(z -I- Agc) Ad 

We have recently argued ^2§\ that theoretical PB distri- 
butions should be convoluted with the experimental res- 
olution function R{z) for comparison with experiments, 
as follows, 

n+{z) ^ J n+{z')R{z' - z)dz'. (6) 
In the particular case of a gaussian resolution function 
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the convolution is 



n+{z) = [ n+{z')e'^-^dz. (7) 
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T he wid th in the Gaussian function is given as a sum, F = 
a/^TI, where each F^'s account for surface-roughness 
due to capillary waves and imperfections in the mono- 
layer and other contributions. Since surface-roughness is 
the dominant contribution we use F w ^Ref, where ^Ref 
is the average roughness determined from the reflectiv- 
ities. From the statistical mechanical perspective, the 
above convolution maybe justified if the fluctuations of 
the interface are gaussian and the wavelength of the cap- 
illary waves A is large (A ^ S, where S is typical ion size) 
so that the ions "see" an effective flat charged interface. 
Alternative methods for computing distributions next to 
flexible interfaces have recently been proposed ,29J . Fig- 
ure ^B) shows convolution of the distributions shown 
in Fig. n^A) for two different values of F as indicated. 
It is interesting to note that the convoluted calculations 
are practically indistinguishable as a function of bulk salt 
concentration. In Section [3 we discuss quantitatively to 
what extent the experimental resolution limits our ability 
to compare with theoretical distributions. 

The system we describe in this paper differs from sim- 
ple PB in that one of the ions species, the hydronium 
H3O"'', can bind to the interface, and as a result, out 
of the Np phosphate groups forming the interface, only 
aNp are electrically charged. The density of charge in 
the system (including the interfacial charge 2; = 0) is 
therefore 



qana{z) - ean(,S{z). 



(8) 



The index a runs over the different species of ions present 
in solution and Qa — ±1 is their valence (i.e. the solu- 
tion contains monovalent ions only). In the free energy, 
the ions are treated as an ideal gas within the solvent, 
and subjected to electrostatic contributions (the stan- 
dard PB) 30]. We also incorporate two additional terms 
that represent the favorable binding of hydroniums to 
the interface and the mixing entropy of the charged and 
neutral phosphate groups at the interface. 



+ kBTj2 [ d'rn,(r)(log(wna(r))-l) 

a •' 

- Eh{1 - a)Np - ksTlog 



(9) 



Npl 



(aiVp)!((l - a)Np)\ J ' 



where Eh is the gain in free energy for a hydronium 
binding to the phosphate group. The quantity v has di- 
mensions of volume and defines the standard state. 

The condition that the free energy is a minimum 
for variations of the different number densities ^^(r) 
under the constraint of a fixed number of particles 



(canonical ensemble) leads to the Boltzmann distribu- 

tion na{z) = v e^s^e ''b^ with chemical potential 
pa = kBT\og{n^v), and also to the Poisson equation Eq. 
(jSJ. Further minimization of the free energy with respect 
to the parameter a yields the additional equation 

EH/kBT-cbiO)+log{aNp)-log{{l-a)Np) = -pn/kBT 

(10) 

where ip{0) is the value of the electric potential at the in- 
terface (the contact value) . This equation may be rewrit- 
ten as 



1 - a 



(11) 
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and expresses the equilibrium of the process HsO^ 
< — > P04H°+Il20 J^. If the standard state is defined 
as V ~ 1 molar, the binding free energy Eh is related to 
the pKa of the molecules forming the interface by Eh = 
knT log Ka and the hydronium concentration is related 
to the pH of the solution by uh = lO^^'^. Substitution 
of these values into Eq. (|ll|l yields 



= IQ-P^" 



1 - a 

The fraction of charged phosphate groups is then 

1 



1 + 10- ipH-pK a) e-'P{o) ■ 



(12) 



(13) 



As a result, the PB equation Eq. (jSJ needs to be solved 
with the self-consistent boundary condition Eq. H13|l 
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1 + io-(pff-p^o)e-0(o)' 



(14) 

The ion distribution is given by Eq. Q with a renormal- 
ized Guoy-Chapman length A'^p — \gcI<^ (RPB). At 
high surface charge, the contact point potential is signifi- 
cantly larger than ksT, and the consequences of Eq. H14() 
are dramatic. The ion distribution becomes strongly de- 
pendent on ion concentration for points very close to the 
interface, as shown in Fig.^| This should be contrasted 
with Fig. ^ where the ion distribution is independent of 
bulk concentration. 



III. EXPERIMENTAL DETAILS 

To manipulate ion bulk-concentrations, we used Csl 
(99.999%, Sigma Corp., Catalog No. 203033) solutions in 
ultra-pure water, taking advantage of the L3 resonance 
of Cs ions at 5.012 keV in anomalous reflectivity mea- 
surements. To control surface charge density, monolay- 
ers of dihexadecyl-hydrogen-phosphate (DHDP, see Fig. 
m (C32H67O4P; MW = 546.86, Sigma Corp., Catalog No. 
D2631) were spread from 3:1 chloroform/methanol solu- 
tions at the air-water interface in a thermostated Lang- 
muir trough . DHDP was chosen for this study, since 
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it forms a sim ple in-p lane structure at high enough sur- 
face pressures [SJ, H and its hydrogen-phosphate head- 
group (i?-P04H) has a pK^ — 2.1, presumably guaran- 
teeing almost complete dissociation [P0j]/[i?-P04H] « 
0.99999, with one electron-charge per molecule (cto ~ 0.4 
C/m^). Monolayers of DHDP were prepared in a ther- 
mostatic, solid Teflon Langmuir trough, and kept un- 
der water-saturated helium environment. Ultrapure wa- 
ter (NANOpure, Barnstead apparatus; resistivity, 18.1 
Mficm) was used for all solution preparations. Mono- 
layer compression, at a rate of ~ 1 A^/ (molecule xmin.), 
was started 10-15 minutes after spreading to allow sol- 
vent evaporation. The monolayer was then compressed 
to the desired surface pressure. During the compression, 
the surface pressure was recorded by a microbalance us- 
ing a filter-paper Wilhelmy plate. The Langmuir trough 
is mounted on a motorized stage that can translate the 
surface laterally with respect to the incident beam to al- 
low the examination of different regions of the sample to 
reproduce results and monitor radiation damage of the 
monolayer. 

X-ray studies of monolayers at gas/ water interfaces 
were conducted on the Ames Laboratory Liquid Surface 
Diffractometer at the Advanced Photon Source (APS), 
beam-line 6ID-B (described elsewhere [s^)- The highly 
monochromatic beam (16.2 and 5.012 keV; A = 0.765334 
and A = 2.47374 A), selected by a downstream Si dou- 
ble crystal monochromator, is deflected onto the liquid 
surface to a desired angle of incidence with respect to 
the liquid surface by a second monochromator [Ge(220) 
and Ge(lll) crystals at 16.2 and 5.012 keV, respectively] 
located on the diffractometer [s^ [s^ . Prior to the mea- 
surements, the absolute scale of the x-ray energy was 
calibrated with six different absorption edges to better 
than ±3 eV. 

X-ray reflectivity and GIXD are commonly used to de- 
termine the monolayer structure [^E 1^ US ■ Spec- 
ular XR experiments yield the electron density profile 
across the interface, and can be related to molecular ar- 
rangements in the film. Herein, the electron density pro- 
file across the interface is extracted by a two-stage refine- 
ment of a parameterized model that best fits the mea- 
sured reflectivity by non-linear least-squares method. A 
generalized density profile p{z) = p'{z) + ip"{z) of the 
electron-density (ED) and the absorption-density (AD) 
(real and imaginary parts, respectively) is constructed 
by a sum of error functions as follows: 

P(^) = 2 ^ '^^^'^ \ ^/^) ^ ^'^^^ ^ ^^'^^ 

where N-l-1 is the number of interfaces, pi = + ipi"\ 
p\ and Pi' are the ED and AD of ith slab, zi and are 
the position and roughness of ith interface, respectively, 
p\ is the electron density of the solution (« 0.334 e/A^ 
), and pAr+2 = is the electron density of the gaseous 
environment. The use of a different roughness for 
each interface preserves the integral of the profile along 



Z or the electron density per unit area, thus conserving 
the chemical content per unit area a. Although small 
variations are expected in f^i for interfaces that separate 
rigid portion of a molecule (hydrocarbon-chains/gas in- 
terface and hydrocarbon-chain/headgroup interface, for 
instance), somewhat larger variations can occur at differ- 
ent interfaces (such as, gas/hydrocarbon chains interface 
versus headgroup/subphase interface). The AD profile is 
particularly important at the Cs resonance (5.012 keV) as 
demonstrated below. The reflectivity is calculated by re- 
cursive dynamical methods [iol I4H of the discretized ED 
and AD in Eq. (|15|) . In the first stage of the refinement, 
the variable parameters used to construct the electron 
density across the interface /9(z) are the thickness values 
of the various slabs di = — Zi|, their correspond- 

ing electron densities pi , and interfacial roughness . By 
non-linear square fit (NLSF) we determine the minimum 
number of slabs or 'boxes' required for obtaining the best 
fit to the measured reflectivity. The minimum number of 
slabs is the one for which the addition of another slab 
does not improve the quality of the fit, i.e., does not im- 
prove x^- In the second stage, we apply space filling and 
volume constraints Js^, 10,23 to calculate pi by assign- 
ing to each slab a different portion of the molecule, the 
ions and water molecules, to a profile that has the same 
number of slabs as obtained in stage one of the analysis. 
As described in detail in section IIVI the model elimi- 
nates the Pi's as parameters, and is self consistent with 
the molecular structure, providing detailed information 
on the constituents of each slab. 

The GIXD measurements are conducted to determine 
the lateral organization in the film. The angle of the 
incident beam with respect to the surface, a, is fixed 
below the incident and reflected critical angle (ac = 
A(/9sro/7r)^/^; tq = 2.82 x 10~^^ cm, where ps is the elec- 
tron density of subphase) for total reflection, while the 
diffracted beam is detected at a finite azimuthal in-plane 
angle, 20, and out-of plane, (3 (the angle of the reflected 
beam with respect to the surface). Rod scans along the 
surface normal at the 2D Bragg reflections are used to 
determine the average ordered chain length and tilt with 
respect to the surface normal. The intensity along the 
rod of the 2D Bragg reflection is analyzed in the frame- 
work of the distorted wave Born approximation (DWBA) 

I{Q.y, Qz) « \t{kz,,)\^\F{Qz)\'\t{kzj)\', (16) 

where t{kz,i) and t{kzj) {kz,i = fcosinQ;;fczj — A:osin/3) 
are the Fresnel transmission functions, which give rise to 
an enhancement at the critical angle. In modeling the rod 
scans, the length and tilt of the tails are varied, exam- 
ining two tilt directions: one toward nearest neighbors 
(NN) and the second toward next NN fNNNl fsa li^ . 
The form factor for the tails is given by 

F(Q:) = sin(QlV2)/(QlV2) (17) 

where is defined along the long axis of the tail, and I 
is the length of the tail. 
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FIG. 2: Surface pressure versus molecular area for DHDP 
spread on Csl solutions at various concentrations as indicated. 
Reflectivity and GIXD were performed at constant surface- 
pressures 30 mN/m and 40 mN/m. 



IV. EXPERIMENTAL RESULTS 



FIG. 3: GIXD scans versus the modulus of the inplane mo- 
mentum transfer Qxy, at surface pressure tt = 30 mN/m 
(curves are shifted by a decade for clarity). The Bragg peaks 
are independent of bulk salt concentration indicating no sig- 
nificant change in in-plane molecular packing. GIXD scan 
for 10"^ M Csl (without DHDP) shows a broad peak at 
QxY ~ 2.0 A~^, due to the surface structure of water. 



A. Isotherm Comparisons 

Surface pressure versus molecular-area (tt — A) 
isotherms of DHDP at various Csl salt concentrations 
(uh) are shown in Fig. |21 For tt > 0, the isotherm 
exhibits two distinct slopes, associated with crystalline 
tilted and untilted acyl-chains with respect to the surface 
normal. The transition from tilted to untilted chains at 
(AtjTTf), occurs at a constant At « 41.5 A^, whereas tt^ 
increases with salt concentration nf,. Csl and other elec- 
trolytes (NaCl and CsCl) in solution significantly influ- 
ence the isotherm causing an increase of the monolayer- 
coalescence area Ac, (i.e., tt > 0) with the increase in 
rife. For A < 39 A^, approximately the cross-section of 
the two acyl-chains of DHDP, (constant tt « 55 mN/m) 
the monolayer is in the yet poorly characterized state of 
collapse. In the present study we focus on the untilted 
crystalline phase (30 ^ tt < 45 mN/m), where the molec- 
ular area variation at a fixed tt is less than 1.5% (in the 
data analysis we allow ± 5% variation in molecular area). 



B. GIXD and Reflectivity off Resonance 

GIXD experiments provided additional insight into the 
molecular packing of the acyl chains within the Lang- 
muir monolayers, namely, the average inplane density of 
the headgroups and the surface charge-density. Figure 
shows GIXD scans from DHDP on pure water and 
on Csl solution (10~^ M, at tt = 30 mN/m) and from 
bare surface of Csl solution (10"'^ M) before spread- 
ing the monolayer as a function of in-plane momentum 



transfer Qxy = \J Q\ + Q\ (ill' The broad peak cen- 
tered at Qxy ~ 2 A"-'^ is due to the structure- factor 
of the aqueous-solution interface. The spreading and 
compression of the monolayer slightly modifies the wa- 
ter structure factor peak at Qxy ~ 2 A~^, and brings 
about two prominent Bragg reflections due the ordering 
of acyl-chains superimposed on a modified surface liquid 
structure- factor 32] . The main features of the diffraction 
pattern are independent of ionic concentration, consis- 
tent with the isotherms at the 30-40 mN/m region, that 
show very small variations in the molecular packing. Fig- 
ure0]shows the 2D diffraction pattern consists of a strong 
Bragg reflection at Qxy — 1.516 A^^ and a weaker peak 
at Qxy = 2.627 A^^, corresponding to 4.145, and 2.392 
A d-spacings, respectively (Table IJ). The shape, spac- 
ing, and location of the two intense peaks correspond 
closely with literature values for (1,0) and (1,1) planes in 
a hexagonal unit cell of ordered alkyl chains, also con- 
flrmed by the ratio i)/Qxy(l, 0) = V^. This 

unit cell (molecular area 19.83 A2) agrees with the cross- 
sectional area of alkyl chain Each headgroup has 
two alkyl chains, giving a molecular area of 39.66 A2, m 
agreement with values obtained from the tt — A isotherm 
40.5 A^, TT = 40 mN/m). The small variations in 
molecular areas are attributed to the existence of domain 
boundaries, defects and minute impurities. The peaks in 
Fig. 01 were fltted to lorentzians the parameters of which 
are given in Tabic The peak linewidth A is signifi- 
cantly larger for the higher order peak. This is expected 
in simple 2D crystals |43 |. and is even more pronounced 
for 2D crystals fluctuating in 3D space (fluctuating teth- 
ered membranes, see Ref. 49]). The inset in Fig. 0| 
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TABLE I: Best-Fit Parameters to high-resolution diffraction 
scan of DHDP on 10"^ M Csl solution (tt = 30 mN/m) 



peak 


QxY (A-1) 


A (A-) 


Intensity (a.u) 


(1,0) 


1.516±0.003 


0.026 


0.522 


(1,1) 


2.627±0.024 


0.060 


0.009 




FIG. 4; Background subtracted GIXD pattern for a DHDP 
monolayer on 10~^ M Csl solution (tt = 30 mN/m) and the 
corresponding rod scan (shown in the inset) at the (1,0) peak 
[QxY = 1.516 A-i). 
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shows the rod-scan at Qxy (1,0) Bragg reflection. The 
the rod-scan analysis (using Eqs. and ((T7I) ) yields 
an average chain-length ^ 20 A, and a tilted angle with 
respect to the surface normal < 5°, consistent with the 
reflectivity and previous reports 34]. 

Figure E^A) shows the normalized reflectivity curves, 
R/Rf (where Rp is the calculated reflectivity of an ide- 
ally flat subphase interface), for DHDP (tt = 40 mN/m) 
on pure H2O and Csl solutions measured at E = 16.2 
keV. The solid lines are the best-fit calculated reflectivi- 
ties based on the ED profiles shown in Figl^J^B). Similar 
reflectivity curves were obtained for tt = 30 mN/m. In 
Fig. [S{A) , all x-ray reflectivity curves differ in the exact 
position and the sharpness of their minima, and the in- 
tensities of their maxima. Given that, as already shown, 
the packing of DHDP is basically independent of salt 
concentration for tt = 40 mN/m, the reflectivity curves 
in FiglSj^A) qualitatively show a strong dependence of 
ion distribution close to the interface on bulk ion concen- 
tration, in quantitative agreement with RPB theory as 
discussed above. 

In the first stage of the analysis, we find that the 3- 
slab {N = 3) model provides good quality flt to all re- 
flectivities and it does not improve with the addition of 
more slabs, i.e., more parameters. Our measured reflec- 
tivity for DHDP on pure H2O is consistent with pre- 
vious measurements ^33| but extends to larger momen- 
tum transfers {Qz), allowing for a more reflned struc- 



FIG. 5: (A) X-ray reflectivity (circles) and corresponding 
best fit (solid lines) for the DHDP monolayer at four solutions 
(tt = 40 mN/m) (curves are shifted by a decade for clarity). 
(B) ED profiles used to calculate the fits shown in (A). 



tural analysis. Indeed, the two-slab model used in Rcf. 
|33| was found to be slightly inadequate, particularly at 
large Qz and a better flt is achieved by adding an ex- 
tra slab at the water-headgroup region, as shown in Fig. 
Thus, our detailed analysis of DHDP on pure wa- 
ter, differs from the one reported in Ref. |33 in which 
the headgroup resides on a thin layer (4 - 6 A thick) of 
ED that is just slightly larger than that of bulk water 
(see Table mi . Similar observations were also reported 
for other monolayers at gas/water interface, and were in- 
terpreted as interfacial water restructuring induced by 
hydrogen bonds. '42^ Further evidence of water restruc- 
turing at the interface is also found in the overall GIXD 
of the interfacial structure-factor of water, especially a 
decrease in peak intensity at Qxy ~ 2 is observed 
[3^ . Modeling DHDP monolayers on the salt solution is 
slightly more complicated as Cs"*" concentration decays 
slowly as a function of distance from the interface. As 
sketched in Fig. [T] we assume that Cs ions are present 
in both the head-group slab and the slab contiguous to 
it toward the bulk of the solution. Table ^ shows the 
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FIG. 6: (A) Reflectivity (circles) taken from DHDP on pure 
H2O and the best fit by using two-slab model (dashed line) 
and three-slab model (solid line). (B) ED profiles extracted 
from two-slab and three-slab model. 
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FIG. 7: Schematic illustration of the three-slab model used 
to calculate self consistently the electron-density profile as- 
suming the a DHDP monolayer of known average molecular 
area, from GIXD and tv — A isotherm, and the associated Cs""" 
and water molecules in the different slabs. Volume constraints 
were also applied in the ED calculations of the different po- 
tions of the molecule and the ion-distribution. 



TABLE II: Best-Fit Parameters to the measured refiectivities 
of DHDP monolayers at tt = 40 mN/m that generate the ED 
profiles across the interface. In this work, the error estimate 
(in parentheses) of a parameter is obtained as described in 
Refs. O 113 by fixing a parameter at different values away 
from its optimum and readjusting all other parameters to a 
new minimum until increases by 50%. Thicknesses of head 
group and Cs box are not well defined due to electron density 
decay from z = to the bulk. 



subphase H2O 



10"^ M Csl 10"^ M Csl 10-^ M Csl 



dtaii (A) 18.7(5) 19.6(10) 19.2(8) 20.2(6) 

ptaii (e/A^) 0.311(8) 0.320(17) 0.329(14) 0.304(8) 

dhead (A) 4.4 3.2 3.8 4.3 

phead (e/A^) 0.476(15) 0.547(27) 0.590(39) 0.624(38) 

dcs (A) 3.6 9.5 3.3 4.2 

PCs (e/A=^) 0.375(11) 0.347(5) 0.431(14) 0.441(8) 



parameters used to produce the ED profiles in Fig.|5jB), 
and the best-fit shown in Fig.E^A). The position at z = 
is defined by the interface between the phosphate head- 
group and the hydrocarbon chain. ED profiles show that 
electron densities at and below the phosphate headgroup 
region are higher with the increase of salt bulk concen- 
tration. The small differences of ED's associated with 
the alkyl-chains for the different subphases are due to 
the minute variations in molecular areas as shown in the 
isotherms above. 



C. Anomalous Reflectivity 

Figure IHl^A) shows refiectivities of DHDP spread on 
10-3 M Csl for TT = 40 niN/m at 16.2 and 5.012 keV. 
The reflectivity taken at 16.2 keV has sharper and deeper 
minima that are slightly shifted to smaller Qz, compared 
to the reflectivity taken at 5.012 keV. This is due to the 
dependence of p'^ and p" on the energy of the x-ray en- 
ergy. At the two energies that the measurements were 
conducted (16.2 and 5.012 keV), p' and p" dramatically 
change only for cesium (as shown in Fig. |3 and only 
slightly for the phosphorous ion, whereas for the remain- 
ing constituents the binding energies are smaller than 
5.012 keV and therefore all electrons behave as free elec- 
trons. Here we apply stage two of the analysis, each 
slab is associated with a portion of the molecule, and the 
ED's and AD's are calculated self consistently applying 
volume constraints. Thus, the ED of the hydrocarbon 
slab is given by 



and 



P'taxl = Ntatl/dtailA 



P'Ll = 



(18) 



(19) 



where Ntaii = 258 is the total number of electrons in the 
two acyl chains [33l|. Similarly, we can calculate pi for 
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FIG. 8: (A) Normalized x-ray reflectivities measured at 16.2 
and 5.012 keV for DHDP monolayer spread on 10^^ M Csl 
solution (tt — 40 mN/m). The solid lines are calculated reflec- 
tivities using the ED profiles shown in (B) . The two data sets 
were combined and refined to a model with common struc- 
tural adjustable parameters. Also shown is the profile of ab- 
sorption factor /3, which at 5.012 keV is dominated by the 
presence of Cs""" close to the interface. (C) Solid smooth line is 
the distribution of Cs""" determined from the refiectivity mea- 
surements. The dotted line is the ion distribution calculated 
from the RPB equation with the corrected surface-charge den- 
sity due to hydronium affinity to PO J . The dashed line is the 
RPB result convoluted with a gaussian of width given by the 
average surface roughness of the monolayer obtained from XR 
without any adjustable parameters. 



the headgroup and for the Cs+ as follows, 

P'l (^Ncs+Zcs+ + NH2OZH2O + ^PO^^PO^^ I'^i^-' 

(20) 

p'l = {pcsNcs+/pOcs+ t^pNpo-/pOp) /2diAXro, (21) 
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FIG. 9: Effective number of electrons and the absorption 
factor f3 for Cs"*". Zeff = p' /n, where n is the number density 
and f3 = \'^p"ro/2n. 



where Nj is the number of ions or molecules, Zj is the 
number of electrons per ion or molecule, fij is the linear 
absorption coefficient of the material when the density 
of material is poj- ^po~~^ headgroup slab and 

N 



po: 



in the Cs slab (see Fig. [Tj). 



Using p-, p", and Eq. l(TB|) . we can create the general- 
ized density p{z) = p'{z) + ip"{z), which includes both 
electron-density and absorption-density. We then apply 
the following volume constraints 



d,A = Ncs+Vc 



Nh^oVh; 



20 



^PO.ypo: 



(22) 



where, Vh^o = 30 A^, y, 



60 A'^ (calculated from 

the reflectivity of DHDP on water), and Vcs+ ~ 20 A^ 
(calculated from the ionic radius in standard tables) ^3 ■ 
The advantage of this method is that a unique set of 
parameters is used to fit both reflectivities at and off 
resonance simultaneously, thus providing a strong self- 
consistency test to the analysis. This is very similar to 
the approach developed to determine the structure of a 
phospholipid monolayer by refining neutron and x-ray 
reflectivities simultaneously [43.l45l |. 

The solid lines in Fig.|SfA) are calculated from the gen- 
eralized density p{z), obtained from parameters of a sin- 
gle model-structure for the combined data sets, as shown 
in Fig.lSjB). The best fit structural parameters obtained 
by this method for various concentrations of Csl in so- 
lution are listed in Table IIIII The absorption factor (3 
shown in Fig. |HlB) can be converted to p", AD curve, 
by a factor (/? = A2p"ro/27r). The AD curve for 5.012 
keV up to a normalization factor is practically the pro- 
file of the counterion Cs+ close to the interface (there is 
a minute contribution to the AD from phosphorous in 
the headgroup region, as shown in Eq. (|21|l 'l. The differ- 
ence between the ED's at and off resonance, normalized 
by Ze//(16.2 keV)-Ze// (5.012 keV),27i], gives the desired 
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FIG. 10: (A) Interfacial Cs^ distributions (solid lines) deter- 
mined Irom anomalous reflectivities (at 16.2 and 5.012 keV) 
for various Csl bulk concentration (shifted by 0.5 M for clar- 
ity). Calculated and convoluted distributions based on RPB 
with renormalized surface charge density as described in the 
text, are shown as dashed lines. (B) Square symbols are num- 
bers of Cs"*" ions per DHDP (« 41 A^) by integrating (up to 15 
A) the experimental distribution obtained from the anoma- 
lous scattering for 10"^, lO"", and 10"^ M. For 10"^ and 
10~^ M, the integrated number of ions are determined from 
the reflectivities off resonance only (the reflectivities at res- 
onance for these concentrations were not measured). The 
dashed line and the solid line are obtained from PB theory 
and RPB theory, respectively. 



ionic distribution at the interface. Figure ISJC) shows 
(solid line) the experimental Cs"*" distribution close to 
the interface at 10~^ M Csl. Similar distributions corre- 
sponding to other bulk Csl concentrations are shown in 
Fig.HniA) (solid lines). 



V. COMPARISON OF EXPERIMENTAL 
RESULTS WITH THEORY 

A. Ion Distributions 



TABLE III: Best-Fit Parameters to the data sets, in which 
the reflectivities measured at and off resonance are combined, 
for various slat concentrations at vr = 40 mN/m. 



subphase(CsI) 


10"^ M 


10""' M 


10"^ M 


10-2 M 


10"^ M 


dtail (A) 


19.6 


19.8 


19.2 


18.8 


20.2 


dhead (A) 


3.2 


3.7 


3.8 


6.7 


4.3 




0.002 


0.013 


0.270 


0.511 


0.523 


dcs (A) 


9.5 


4.7 


3.3 


6.9 


4.2 




0.119 


0.288 


0.289 


0.187 


0.410 


Area (A^) 


41.00 


41.04 


40.97 


41.00 


42.08 


total iVcs+ " 


0.12(5) 


0.30(7) 


0.56(10) 0.70(12) 0.93(12) 


" Number of C 


s"*" in the headgroup slab. 


' Number of Cs+ 



in the Cs slab. estimated errors are given in parentheses. 



experimental distribution along the z-axis and can be 
checked self-consistently from the model used in the anal- 
ysis of the combined data set (Table llTT|) . The number of 
ions per DHDP versus Csl bulk concentration are plot- 
ted in Fig trUT B) with square symbols. In Fig. [TUT B) the 
integrated values obtained by PB theory with the sur- 
face charge corresponding to the fully deprotonated phos- 
phate groups is also plotted (triangles connected with a 
dashed line) for comparison. As shown, the experimental 
integrated number of Cs^ at the interface varies roughly 
as a power-law of bulk concentration, which is well de- 
scribed with RPB (solid line) without fitting parameters, 
as described in Section ITTl 

The ion distribution predicted from RPB theory us- 
ing Eq. with the renormalized Gouy-Chapman length 
(dashed line) is compared with the experimental distri- 
bution in Figs. IHIC) and HHIA) (solid line). As dis- 
cussed in Section |n] the theoretical distribution needs to 
be convoluted with the effective experimental resolution 
function. The distribution resulting from the convoluted 
RPB, with no fitting parameters, with T (F « 3.8 A) 
obtained from the analysis of the reflectivity, is shown 
to reproduce the experimental data remarkably well, as 
shown in Fig. |HIC). The value for pH — pKa = 4.5, used 
in the calculation, is consistent within the range of our 
measured values for the pH ~ 6.5. 

The distributions corresponding to the three bulk Cs"*" 
concentrations 10"^, 10"^, and 10"^ M are shown in Fig. 
llOf A) with a solid line. The agreement with the RPB 
(dashed lines) convoluted as described is remarkable ex- 
cept for points far from the interface. We attribute this 
error to the difficulty to include slowly decaying tails of 
the PB theory to the ED profile modeled by the faster 
decaying Error functions. 



B. Possible deviations from PB theory 
distributions 



We first compute the integrated number of Cs+ per 
DHDP over the first 15 A next to the charged inter- 
face. This number can be obtained by integrating the 



This Section deals with the examination of the sen- 
sitivity of the anomalous reflectivity in determining the 
ion distribution, in addition to the integrated number of 
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FIG. 11: (A) Solid line corresponds to the RPB profile at bulk 
concentration 10~^ M. Step-like functions preserve the inte- 
gral (up the first 15 A) of the continuous RPB distribution. 
(B) Corresponding convoluted distributions with a gaussian 
of width r = 3.8 A. 



ions at the interface. The results presented in Fig. llor A) 
show the unique capability of the anomalous reflectiv- 
ity technique in providing ion distributions. Whereas 
the integrated number of ions, such as the ones shown 
in Fig. llOf B') , can be obtained from standard reflectiv- 
ity measurements and other experimental techniques, as 
discussed in Section ^ the determination of ion distribu- 
tions requires the use of anomalous reflectivity. In fact, 
data such as shown in Fig. llor B) has been used to assess 
the validity of PB theory in the past. As discussed in the 
Introduction, the excellent agreement maybe somewhat 
deceptive in that it may hide significant short-distance 
deviations from RPB theory because only the total inte- 
grated ion density is involved. Our experimental resolu- 
tion constrain possible deviations from RPB distributions 
to within 3.8 A. If such deviations are non-monotonic, 
showing bumps or oscillations, then the constraints are 
even more stringent. As an example, in Fig. Illf A') we 
construct two step- wise distributions (TV = 6, and = 8 
steps of 2.5 and 1.9 A width, respectively) whose to- 



tal integrated area is the same as for RPB (continuous 
line). Although the histograms have the same area as 
the RPB, they incorporate hypothetical non-monotonic 
decay of the distribution in the form of bumps. As it is 
shown in Fig. lllf B). despite the fact that the step-size are 
smaller than the 3.8 A resolution, such distributions can 
be ruled out by the experiments. This demonstrates that 
if actual ion distributions are non-monotonic, their max- 
ima or minima must be short-ranged in nature (shorter 
than « 2 A) to be consistent with our data. 



VI. CONCLUSIONS 

The goal of this study was to explore the accuracy of 
PB theory for monovalent ions. For that, we selected 
a system with a relatively high surface charge density 
(one electron per 40 A^, lattice constant ^ 6.8 A) and 
a low pKa = 2.1 system, which we expected would pro- 
vide the most favorable scenario to observe deviations 
from PB theory. From the results obtained in our exper- 
iments we conclude that PB with the renormalization of 
surface charge density due to proton-transfer and release 
processes (RPB) is strikingly accurate. 

Certainly, the accuracy of our results is limited by 
the effective experimental resolution, which is dominated 
by the natural surface-roughness of the gas-amphiphile- 
water interface (« 3.8 A), due mainly to thermally ac- 
tivated capillary- waves. We should point out, however, 
that this resolution is comparable to the diameter of one 
Cesium ion (3.2 A) or a single water molecule {'-^ 3 
A), and therefore our experimental distributions con- 
straint deviations from RPB theory to very short-range 
variations involving one, at most two, water molecules 
or cesium ions. If the actual ionic distribution is non- 
monotonic, departures from RPB theory are even more 
constrained as it follows from the discussion in Fig. 1111 

Our results show that theoretical effects that are usu- 
ally suggested to modify PB theory, such as finite ion- 
size, in-plane modulations of surface charge density, hy- 
dration forces, short-range interactions, the roughness of 
the surface or image charges are not necessary to describe 
the experimental data. This is not to be understood as 
implying that such effects are not present, but rather 
that their significance is entirely limited to a character- 
istic spatial distance of the order of « 4A or less. As 
for claims based on the modification of PB by hydration 
forces (see for example, Ref. ^2)' ou'" experimental re- 
sults conclusively rule out the possibility of modifications 
of PB within the 10 - 20 A range from the interface and 
confine such corrections, if present, to within the first 
3 - 4 A from the interface as discussed. Although the 
experimental counter-ion distributions are well described 
within the RPB theory, we point out that the reflectivity 
and GIXD hint at water restructuring at the interface. 
Future theoretical or numerical work may clarify this is- 
sue. 

A recent report analyzing the accuracy of PB near 
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charged liquid-liquid interfaces by the use of x-ray sur- 
face sensitive techniques shows that ion distributions are 
well described by PB at concentrations of 10 mM, but 
marked deviations are found at higher concentrations (of 
the order 100 mM), where ion-ion correlations have to 
be included in the theoretical analysis. [H^ We did not 
perform anomalous reflectivity for concentrations above 
lOmM, so we are not able to provide ionic distributions 
for these concentrations. We point out, however, the 
good agreement found for the integrated quantities at 
these concentrations fFig. llur B)). which provides an ex- 
ample of integrated quantities possibly hiding deviations 
from actual distributions as pointed out in the introduc- 
tion. 

The results presented in this study enhance our under- 
standing of the electrostatics in aqueous media, and also 
show the strength of surface sensitive x-ray synchrotron 



techniques in obtaining high resolution data. 
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